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The effects of carrier-impurity correlations due to a Kondo-like spin-spin interaction in diluted 
magnetic semiconductors are investigated. These correlations are not only responsible for a transfer 
of spins between the carriers and the impurities, but also produce non-perturbative effects in the 
spin dynamics such as renormalization of the precession frequency of the carrier spins, which can 
reach values of several percent in CdMnTe quantum wells. In two-dimensional systems, the preces¬ 
sion frequency renormalization for a single electron spin with defined wave vector shows logarithmic 
divergences similar to those also known from the Kondo problem in metals. For smooth electron 
distributions, however, the divergences disappear due to the integrability of the logarithm. A pos¬ 
sible dephasing mechanism caused by the wave-vector dependence of the electron spin precession 
frequencies is found to be of minor importance compared to the spin transfer from the carrier to 
the impurity system. In the Markov limit of the theory, a quasi-equilibrium expression for the 
carrier-impurity correlation energy can be deduced indicating the formation of strongly correlated 
carrier-impurity states for temperatures in the mK range. 
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I. INTRODUCTION 

A perturbative treatment of the interaction between 
quasi-free electrons in a metal with localized mag¬ 
netic impurities predicts logarithmic divergences in sev¬ 
eral quantities such as resistivity and entropy at zero 
temperature^’^. This finding, the Kondo effect, is a fa¬ 
mous example of a situation where perturbation the¬ 
ory leads to unphysical conclusions whereas the mea¬ 
sured values of the resistivity assume finite values. The 
Kondo problem, i. e., the question of how to properly 
describe the low-temperature limit of a system with a 
spin dependent carrier-impurity interaction theoretically, 
has opened up a wide field of phyiscs. Although the 
Kondo problem, as it was originally formulated, has been 
solved^, the Kondo physics experienced a revival since 
it has become possible to study experimentally similar 
problems in other systems, e. g., structures where quan¬ 
tum dots play the role of the magnetic impurities^^^^. 
The common feature in these systems is that a micro¬ 
scopic exchange coupling gives rise to an effective Kondo- 
Hamiltonian that assumes the form of a spin-spin contact 
interaction between the quasi-free carriers and the local¬ 
ized magnetic impurities, or quantum dots, respectively. 

Other systems which are usually modelled by a Kondo- 
like Hamiltonian are diluted magnetic semiconductors 
(DMS) where typically II-VI or III-V semiconductors are 
doped with magnetic impurities, usually Mn which ef¬ 
fectively forms a spin-| system. These materials have 
been studied extensively in the last decades^^^^^ due to 
their optical and magnetic properities which make them 
promising candidates for future spintronics devices^"^^^®. 
While the carrier-impurity correlations play a crucial 
role in the metallic Kondo problem, they are often ne¬ 
glected in studies about DMS by employing a mean-field 
approximation^^. In some articles®^ it is argued that the 


situation in DMS is different from the original Kondo sit¬ 
uation in that in the latter only a few magnetic impurites 
and a huge number of quasi-free carriers are present in 
the metal, whereas in the former case, in particular in the 
case of (intrinsic) II-VI DMS, the number of impurities 
usually exceeds the number of carriers. 

On the other hand, a third-order many-body pertur¬ 
bation theory based on the pseudofermion formalism®® 
reveals Kondo-like divergences in the propagator for the 
spin dynamics in DMS due to the hole-impurity exchange 
interaction. From this it was concluded that the carrier- 
impurity correlations should in fact be important for the 
dynamics in DMS. Furthermore, some studies®®^'*'^ sug¬ 
gest that bound magnetic polarons can play a key role 
for the ferromagnetic behaviour of some DMS. These po¬ 
larons consist of carriers which are bound to the mag¬ 
netic impurities and, thus, the carriers and impurities 
are strongly correlated. 

In this article, we address the question of the impor¬ 
tance of carrier-impurity correlations in DMS and the 
relevance of possible Kondo-like divergences. We base 
our study on a microscopic quantum kinetic theory de¬ 
rived by a correlation expansion scheme^® that is ca¬ 
pable of a non-perturbative description of highly non¬ 
equilibrium situations. One aspect of the effects of the 
carrier-impurity correlations on the spin dynamics has 
already been found in previous works'^®^®®: The corre¬ 
lations mediate the transfer of spins between the carri¬ 
ers and the impurities. Since in the Markovian limit, 
the quantum kinetic theory contains the special case of 
rate equations which can also be derived by a Fermi’s 
golden rule approach^®, this spin transfer can, in fact, be 
treated perturbatively®®. Note that in some situations, 
e. g., for excitations close to the band edge in two- and 
lower-dimensional DMS®^, the Markov limit is not a good 
approximation so that deviations from a golden-rule-like 
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exponential decay are predicted. 

In the present study, we show that the carrier-impurity 
correlations are also responsible for another effect in the 
spin dynamics that is not predicted by a perturbative 
method: a renormalization of the precession frequency 
of carrier spins compared with its mean-field value. It is 
shown that the frequency renormalization also contains 
Kondo-like logarithmic divergences in the Markov limit 
in two-dimensional systems. However, these divergences 
never lead to unphysical results in the spin dynamics. 
This is, first of all, due to the fact that the singularities 
are integrable and yield finite values for a non-singular 
spectral electron distribution. Moreover, the divergence 
in the frequency renormalization is only found for t —>■ oo 
where the amplitude of the precessing electron spin has 
already decayed to zero. The Markov limit of the quan¬ 
tum kinetic theory also allows to find an expression for 
the carrier-impurity correlation energy which shows a 
similar behaviour as the frequency renormalization, in¬ 
cluding Kondo-like logarithmic divergences in the two- 
dimensional case. 

The article is structured as follows: First, the quan¬ 
tum kinetic theory is briefly reviewed as well as ef¬ 
fective (PESO, precession of electrons and correlations) 
equations^® based on the quantum kinetic theory. Then, 
the frequency renormalization described by the PESC 
equations is calculated and compared with the result of 
a Markovian approximation to the PESC equations in 
two and three dimensions. A possible electron spin de¬ 
phasing mechanism due to the wave vector dependence 
of the frequency renormalization is discussed. Finally, we 
investigate the mean carrier-impurity correlation energy. 

II. THEORY 
A. System 


The Hamiltonian for conduction band electrons in 
DMS is modelled by 
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where Hq describes the band structure and Hsd is the 
Kondo Hamiltonian which originates from the exchange 
interaction between the s-type conduction band electrons 
and the d-electrons of the magnetic ions. Throughout 
this article, we assume a parabolic band structure with 
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Wk = ■''^here m* is the effective mass. Jsd and V 
are the coupling constant and volume of the DMS, 
and Co-k are the creation and annihilation operators for 
electrons with spin index cr and wave vector k. R/ is 


the position of the /-th magnetic impurity and = 
\I, n){I, n'\ are the projection operators corresponding to 
the spin state of the /-th impurity, e. g., for spin | Mn 
impurities, n = {—|,—1...|}. S„„/ and Sco-' are the 
spin matrices for spin | and ^ systems, respectively. 

B. Equations of motion 

A microscopic quantum kinetic theory based on a cor¬ 
relation expansion scheme was constructed in Ref. 45, 
where equations of motion have been derived for the elec¬ 
tron and impurity density matrices and as well 
as their correlations which are defined by 
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where the brackets denote the quantum mechanical aver¬ 
age as well as an average over homogeneously distributed 
impurities. The equations of motion for these dynamical 
variables are given in Ref. 48. 

The full quantum kinetic equations are lengthy and 
their solution requires considerable numerical effort. 
However, it was found in Ref. 49 that they can be dras¬ 
tically simplified in the case where the number of im¬ 
purity ions Nmu is much larger than the number of the 
quasi-free electrons Nf,. This is usually fulfilled espe¬ 
cially in H-VI DMS where the magnetic doping with Mn 
does not simultaneously lead to p- or n-doping and the 
carriers stem exclusively from optical excitation. To un¬ 
derstand the effective equations derived in Ref. 49 it is 
instructive to first consider the mean-field dynamics for 
the spin Sk = of electrons with wave 

vector k and the impurities (S) = In 

the mean-field approximation, i. e., if the correlations are 
neglected, one finds 



where ojm ■= Eq- (3b) follows from the 

total spin conservation of the Kondo Hamiltonian. In 
the case Nmu Ne, the change of the impurity spin is 
marginal and can therefore be neglected. The precession 
of the electron spin around the mean field due to the 
impurity magnetization, on the other hand, is in general 
important. Eq. (3a) is solved by 

Sk = /?<S>(wM<)Sk> (4) 

where Rn{a) is the matrix describing a rotation around 
the vector n with angle a and the precession frequency 
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wm = ■ (S)/|(S)| is defined so that it has the same 

sign as the coupling constant Jgd- In the mean field ap¬ 
proximation is constant. However, if we also account 
for the carrier-impurity correlations, changes slowly 
with time and constitutes the electron spin in a rotat¬ 
ing frame. If the correlations are formally integrated and 
inserted into the corresponding equations of motion for 
the electron variables, the effective equations for the elec¬ 
tron spin component perpendicular to the impurity 
magnetization can be given as^®: 




dt 
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The coefficients in Eq. (5) are given by := ± 

6° := and := where the compo¬ 

nent of the impurity spin operator in the direction of 
the mean impurity spin is S'11 := S • and the rele¬ 

vant second moments of the impurity spin operator can 
be separated into parallel (S^l ) and perpendicular parts 
(S-*-^) = 5 (S^ — Sll^). The memory function 
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has to be interpreted as an integral operator and the 
time-dependent variables that appear after in 

Eq. (5) are evaluated at t'. Einally, are the occupa¬ 
tion numbers of the states with wave-vector k, i. e., the 
diagonal elements of the density matrix with respect to 
the spin indices. Eq. (5) together with the correspond- 

•t- / I 

ing equations for given in Ref. 49 are called pre¬ 
cession of electron spins and correlations (PESO) equa¬ 
tions, since besides the electron spin, also the correla¬ 
tions := X]cti<T 2 ■ ^nin 2 Qa-lnlki exhibit 

a precession-like movement around the mean field due to 
the impurity magnetization. Note that Eq. (5) is equiva¬ 
lent to the full quantum kinetic theory of Ref. 45 except 
that some source terms for the correlations are neglected 
that are numerically insignificant [cf. Ref. 49 for details]. 

Eq. (5) is only complicated and numerically challeng¬ 
ing due to the time integral induced by the memory func¬ 
tion . Now, working in the rotating frame allows us 
to assume that the electron variables change only slowly 
in time and can equally well be evaluated at t instead of 


t'. The memory integral consists then only of 
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where V is the Cauchy principal value. The Markov ap¬ 
proximation (7) was established by letting t —>■ oo in 
the lower limit of the integral and using the Sokhotski- 
Plemelj theorem. The validity of the Markovian approx¬ 
imation can in general depend on the values of k, ki, 
t as well as the timescale of the change of the electron 
variables and therefore has to be checked numerically. 

If only the real part of the memory function is used 
in Markov approximation and the imaginary part is ne¬ 
glected, the PESC-equations assume a golden rule-type 
form, where the spin transfer dynamics follows approx¬ 
imately an exponential decay to the equilibrium value 
with rate 


E) 


-1 


hW 

-\- D{u!i)b'^'^ 


^b+ ^b- 

D[wi — ujm)-^ + D[wi -\- 

( 8 ) 


for an electron with kinetic energy wi, if the terms of 
second order of the electron variables in Eq. (5) are 
neglected^®. In the expression for the rate, D{(jj) de¬ 
scribes the spectral density of states and depends on the 
dimensionality of the system. 


C. Frequency renormalization in the Markov limit 

One issue that we would like to focus on in the present 
work is the change in the precession frequency described 
in Eq. (5) by the terms proportional to the imaginary 
part of the memory function. Such a renormalization 
of the precession frequency would be absent in any trun¬ 
cated perturbative approach^^. It originates, like the spin 
transfer described by the real part of the memory func¬ 
tion, from the carrier-impurity correlations. 

It is noteworthy that the frequency renormalization is 
singular in the Markov limit described in Eq. (7), i. e., 
the imaginary part of the memory function Gu>^f diverges 
if Wk = Wki- However, this divergence does not lead to 
an unphysical behaviour. Eirst of all, the divergence is 
a feature of the Markovian limit. For finite times t, the 
1. h. s. of Eq. (7) is a finite integral over an analytic 
function and is therefore also analytic. For Wk = Wki, 
the value of the integral is t which only goes to infin¬ 
ity in the Markov limit. As only the electron spin com¬ 
ponent perpendicular to the impurity magnetization is 
affected by the frequency renormalization and this com¬ 
ponent decays approximately exponentially to zero, an 
infinite precession frequency is never observable. 

Similar to the Markovian spin transfer rate in Eq. (8), 
an expression for the frequency renormalization Aw can 
be given in the Markov limit of Eq. (5), if the imaginary 
part of Eq. (7) is used: 
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where, for the sake of simplicity, the terms propor¬ 
tional to in Eq. (5) were neglected, since they 

only matter if a large number of carriers is present. In 
two-dimensional systems, the spectral density of states 
D'^'^{uj) = i® constant, where A is the sample 

area and 0(x) is the step function. In three dimensions, 

a/w ©(w) is proportional to the 

square-root of w. The corresponding frequency renormal¬ 
izations are: 
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where d = Vj A is the quantum well width, and 
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It should be noted that in two and three dimensions 
the frequency renormalization depends explicitly on the 
frequency ojbz, which corresponds to the energy at the 
end of the first Brillouin zone, and diverges in the limit 
ujbz —>■ oo. For typical pump-probe experiments with 
diluted magnetic semiconductors, carriers are optically 
excited relatively close to the band edge. For the excited 
electrons, one can safely assume wi ± ojm ‘C ujbz- In 
this case, we find a logarithmic dependence on ujbz in 
the two-dimensional frequency renormalization. 

With the same assumption also the integral in 
Eq. (10c) for the three-dimensional renormalization can 
be simplified to 


duj- 


w — Wo 


2.yj0JBZ — 7’‘\/|wo|0(—Wo). (11) 


Thus, we find a square-root dependence of the frequency 
renormalization on the cut-off frequency ojbz, as well as 
a square-root dependence on wo = wi ± ojm which only 
contributes if wo is negative. 

The divergence in the limit ojbz —>■ oo is similar to the 
metallic Kondo effect where the divergence in the resis¬ 
tivity is also logarithmic in the bandwidth^. The Kondo 
problem in metals resembles rather the two-dimensional 
than the three-dimensional case in DMS, because the for¬ 
malism for the solution of the Kondo problem usually de¬ 
scribes the carrier system as possessing a constant spec¬ 
tral density of states^. 

It is noteworthy that the divergence of the frequency 
renormalization for wq = wi ± ojm vanishes in the three- 
dimensional case due to the integral over the density 
of states. In two-dimensional systems, a diverging fre¬ 
quency renormalization remains, but only for electrons 
with a unique value of the kinetic energy. For realistic 
optical excitation, however, a smooth spectral electron 
distribution can be expected so that the change of the 
total precession frequency comprises an averaging over 
frequency renormalizations of nearby states. Since the 
logarithmic divergence is integrable, the total frequency 
renormalization remains finite. 


III. NUMERICAL CALCULATIONS 

In order to check the validity of the Markov approxima¬ 
tion for the renormalization of the precession frequency 
of the electrons, we compare the Markov result with cal¬ 
culations, where the memory is taken into account ex¬ 
plicitly. It seems straightforward to use Eq. (5) with 
the time-integral operator Gj^ defined in Eq. (6) and 
solve the integro-differential equations numerically. This 
is, however, a very challenging problem for the following 
reasons. 

From the Markovian expression for the frequency 
renormalization, we find the explicit dependence on the 
value of the cut-off energy Hojbz- Therefore, also oscilla¬ 
tions with frequencies close to ujbz have to be resolved, 
which are on the timescale of a few fs since hwBZ is in the 
eV range. On the other hand, relevant changes of the to¬ 
tal electron spin takes place in the 10-100 ps range. Fur¬ 
thermore, for each time step the calculation of each 
requires a sum over all possible k-states so that the prob¬ 
lem has the complexity 0{N‘^) where Nk is the number 
of discretization points for the k-space. Note that also 
in /c-space, the details of excitations close to the band 
edge in the meV range as well as the full Brillouin zone 
up to energies of a few eV have to be resolved. Such a 
problem also arises in the metallic Kondo effect where 
numerical procedures, such as the famous renormaliza¬ 
tion group®^, have been developed to deal with the large 
value of the band width^. Note that solving the integro- 
differential equation by finding an auxiliary variable, so 
that the problem can be transformed into an ordinary 
differential equation, is equivalent to using the original 
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quantum kinetic theory^®. 

Here, we solve this problem by using approximations 
that allow a separation of electron spins with different 
wave vectors, so that we find a 0{Nk) problem for an 
individual electron with wave vector ki. First of all, it 
is noteworthy that in Eq. (5) couples to the occupa¬ 
tions of states with different wave vectors k. These 
terms, however, are of second order in electron variables 
and have a marginal effect on the dynamics of the per¬ 
pendicular spin component^®, especially if the electron 
density is small, as is usually the case for optically ex¬ 
cited carriers. Neglecting these terms, we can formulate 
equations of motion for the complex perpendicular elec¬ 
tron spin variable (in the rotating frame): 
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where it is assumed that the impurity magnetization 
points in the z-direction. Then, the PESC-equations (5) 
assume the form: 
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It can be seen immediately from Eq. (13) that in the 
equation for electron variables of states with other 
wave vectors only enter in the last term, i. e., the term 
proportional to sj^(t') -I- (t'). Note that a time integra¬ 
tion of cos[(a;k — Wki)(t' — t)] yields which 

has a pronounced peak at Wk = Wki ■ Thus, if the electron 
spin distribution is assumed to be a smooth function in k- 
space, the main contribution of the last term in Eq. (13) 
will be approximately the same if we set 
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u, ={1,-1,1,-1}, (15c) 

X,= (1,-1,0,0}, (15d) 

6^ = (6+y-, fell, fell} (15e) 


and initial conditions = 0 for t = 0. Calculating 

the dynamics of a single electron spin using Eqs. (15) has 
the complexity 0{Nk) and can be done without the need 
for a numerical renormalization group procedure. 


IV. RESULTS FOR THE FREQUENCY 
RENORMALIZATION 

The parameters used for the numerical calculations de¬ 
scribe a Cdo.gsMno.orTe sample with coupling constant 
Jsd = —15 meVnm^, effective mass m* = 0.093 mg®®, 
where mo is the free electron mass, and, in the case of a 
two-dimensional system, a quantum well width of d = 5 
nm. The cut-off energy was taken to be Hlobz = 3 eV. 
The initial impurity magnetization was modelled to be 
thermally distributed and is therefore completely defined 
by the mean value (S'!!) G [—§;§]■ 

We assume that electrons have been spin selectively 
prepared by optical excitation with circularly polarized 
light so that the initial electron spin is perpendicular 
to the initial impurity magnetization (Voigt geometry). 
Eqs. (15) are used to calculate the finite-memory spin 
dynamics for electrons with a defined wave vector ki or, 

equivalently, kinetic energy Huji = 2 ;;^. An exponen¬ 
tially decaying cosine 

sZ (i) ~ (O)e-G^^ cos(uj'j^t) (16) 


4(0 « 4 ( 0 - ( 14 ) 

This approximation was shown to reproduce the non- 
Markovian features of the spin transfer in Ref. 52. Also, 
in contrast to the other terms, the last term of Eq. (13), 
where the approximation is used, does not influence the 
frequency renormalization, due to the absence of an imag¬ 
inary part of the oscillating prefactor cos[(a;k — — 

t)]. Now, with the help of approximation (14), we end up 
with completely decoupled equations for the spins of 
electrons with different wave vectors ki. 

Finally, it is useful for the numerical solution of the 
integro-differential equation (13) to transform it into an 
ordinary differential equation using auxiliary variables 


is fit to the non-Markovian spin dynamics in order to 
find a value for the effective decay rate tJ^(wi) and the 
precession frequency The relative renormaliza¬ 

tion of the precession frequency is given by with 
Aw = uj'jyj - UJM- 

Fig. 1 shows the relative frequency renormalization ob¬ 
tained from a fit to the non-Markovian calculation and 
the corresponding Markovian result for a d-like initial 
spectral electron distribution as a function of the kinetic 
energy hiui. First of all, it can be seen that in three- 
dimensional as well as in two-dimensional systems the 
Markovian and non-Markovian results coincide. In the 
three-dimensional case, the square-root energy depen¬ 
dence of the renormalization for wi < com can be seen 
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FIG. 1. Dependence of the relative frequency renormaliza¬ 
tion on the kinetic electron energy for two- and three- 
dimensional systems according to the calculation including a 
finite memory [Eqs. (15)] (points) and in the Markov limit 
[Eqs. (10)] (lines) for different values of the average impurity 
spin ](S)] (in units of ft). 
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FIG. 2. Spin transfer rate in Markov approximation (line) 
according to Eq. (8) and exponential fit to the calculation 
including a finite memory (crosses) using Eq. (15) for ](S)] = 
0.05 in a two-dimensional system. 


clearly, while in two dimensions, the logarithmic diver¬ 
gence at LOi = ujm is apparent. The positive relative fre¬ 
quency renormalization in three dimensions describes an 
increase in the modulus of the precession frequency. In 
two dimensions, the slightly positive background of the 
renormalization is overcompensated by a negative value 
in the region around the divergence. 

In Fig. 2, the spin transfer rate according to the 
Markov approximation is compared with the value ob¬ 
tained by the exponential fit to the non-Markovian result 
for a calculation with | (S) | =0.05 in two dimensions. The 
step in the rate at wi = wm, which is predicted in the 
Markov limit [cf. Eq. (8)], is found to be slightly rounded 
off in the non-Markovian calculation, but the deviations 


FIG. 3. Relative frequency renormalization in a two- 
dimensional system for a Gaussian spectral electron distribu¬ 
tion centered at Ec = huiM with standard deviations of Es = 1 
meV and Es = 0.1 meV, respectively. The initial electron dis¬ 
tribution as a function of the kinetic energy is visualized in 
the inset as the blue dash-dotted line {Es = 1 meV) and green 
dotted line {Es =0.1 meV) together with the corresponding 
frequency renormalization for d-like excitations (red line) for 
](S)1=0.05. 


between both results are rather small. 

In order to find an estimate for the strength of the 
change of the precession frequency for a more realistic 
electron distribution. Fig. 3 shows the relative precession 
frequency renormalization as a function of the average 
impurity spin where the initial spectral electron distri¬ 
bution [cf. inset of Fig. 3] was assumed to be Gaussian 
with center at Ec = JEjm and standard deviation Eg = 1 
meV (0.1 meV) corresponding to a full width at half max¬ 
imum (FWHM) of ~ 2.35 meV (0.235 meV) or a Gaus¬ 
sian envelope of an exciting laser pulse with a duration 
(FWHM) of ~ 140 fs (1.4 ps). The calculations for Fig. 3 
were performed using the 2D Markovian expression for 
the rates in Eq. (8) and the renormalized precession fre¬ 
quencies in Eq. (10a). It can be seen that the magnitude 
of the frequency renormalization can reach values of sev¬ 
eral percent of the mean-field precession frequency and 
is negative for small values of |(S)|. For larger values of 
the impurity magnetization, the frequency renormaliza¬ 
tion approaches a small positive value. One could expect 
that the narrower electron distribution {Eg =0.1 meV) 
is closer to the d-like case than the wider distribution 
{Eg = 1 meV) and therefore the frequency renormaliza¬ 
tion should be more pronounced. However, it can be 
seen from Fig. 3 that this is only the case for very low 
values of |(S)| (below 0.01 in the case studied here). For 
higher values of the impurity magnetization, the relative 
frequency renormalization approaches the positive back- 
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FIG. 4. Rate obtained from an exponential fit for Gaus¬ 
sian initial spectral electron distributions (cf. Fig. 3) as a 
function of the average impurity spin with (blue dash-dotted 
line/green dotted line) and without (red solid line) accounting 
for a renormalization of the precession frequency. 

ground much faster in the calculations with the narrower 
electron distribution. 

Note that in order to be able to measure or fit a pre¬ 
cession frequency, at least one period of the oscillations 
should be visible before the spin polarization is decayed. 
Thus, the minimal value of the impurity magnetization, 
where one can reasonable deduce a precession frequency 
from the time evolution of the spin polarization, is given 
by \wm\ ^ which yields, for the parameters above, 
|(S)| > 0.01. Therefore, we find that short laser pulses 
with pulse durations of the order of 100 fs provide the 
most promising configuration for experiments to measure 
the frequency renormalization. 

Since the frequency renormalization depends on the ki¬ 
netic energy and therefore the wave-vector of an electron, 
the question arises, whether this dependence leads to a 
dephasing of spins of electrons with different k-vectors. 
To address this question, we show in Fig. 4 the value of 
the rate obtained by an exponential fit to the time 
evolution of the total carrier spin polarization, where the 
same Gaussian initial electron distributions are used as 
in Fig. 3. It can be seen that calculations, where the cor¬ 
relation induced frequency renormalization is neglected. 


produce very similar decay rates as calculations that ac¬ 
count for this renormalization for most of the possible 
values of the impurity magnetization. Only in a regime 
where the impurity spin is small we find a slightly larger 
value (< 1%) of the rate at |(S)| a; 0.05 for Eg = 1 meV 
and |(S)| Ri 0.005 for Eg = 0.1 meV. This increasing de¬ 
cay is the consequence of the dephasing of electron spins 
due to the k-dependence of the frequency renormaliza¬ 
tion. Since the expression for the rate in the Markov 
limit [cf. Eq. (8)] and the frequency renormalization [cf. 
Eq. (10a)] depend on the same parameters, this dephas¬ 
ing mechanism is always accompanied by a genuine spin 
transfer between impurities and carriers that is typically 
much faster than the dephasing itself. 

V. CORRELATION ENERGY 

Most studies on DMS which probe the energies of elec¬ 
trons in DMS use the mean-field approximation^^ and 
describe the effects of the impurity magnetization as a 
renormalization of the electron g-factor which is known 
as the giant Zeeman effect^®. If, however, the build-up 
of carrier-impurity correlations is taken into account, the 
mean s-d exchange interaction energy {Hsd) will devi¬ 
ate from the mean-field value. The correlation energy 
can, in principle, have an impact on the thermodynamic 
properties of DMS which could help, e. g., in the descrip¬ 
tion of the paramagnetic-ferromagnetic phase transition 
in GaMnAs. 

Since the derivation of the PESC-equations in Ref. 49 
required finding explicit expressions for the correlations, 
we can use this theory to get the correlation induced 
correction to the mean-field exchange interaction 

energy analytically: 

{Hsd) = ^ E = 

Inn' 

crcr^kk^ 

=: E • Sk + (17) 

k 

Using the time-integral form of the correlations from 
Ref. 49 we find: 


i^sd ) — ly' ^Mn EE^Sk =-^E^ 


kk' Q!=l 


kiko 


+ 




^ - Y’^ki + + ’^L”ki)] +4^) - +"^^ 1 )) 


(18) 


To understand Eq. (18) it is important to recall that the correlations typically build up on the timescale of a few 













, while the spin-up and spin-down occupations change 
on a ps timescale'^®’'*®. Thus, can be interpreted 

as a quasi-equilibrium value of the correlation energy for 

■f /I 

given values of adiabatically changing occupations . 

As in the discussion of the frequency renormalization, 
we neglect terms of second order in the electron variables 
and apply the Markov approximation to find for the two- 
dimensional case: 


(ffCor\ ^ _ ^sd ''^Mn Am* ^ 

^ - fi V 2Trh ^ 


In 
-fin 
+ ln 


WbZ — (Wki + Wm) 


Wki -f Wm 
WBZ — (wki — Wm) 


b nl^ + 


b'^ni 


Wki — Wm 
CObZ — Wki 


Wki 


y(4i +«ti) 


(19) 


The mathematical structure of the correlation energy 
i^siT) ^9- (19) is very similar to that of the frequency 
renormalization in Eq. (9). To see this relation, it is help- 

ful to express the occupations of the spin-up and 
spin-down subbands in terms of the occupation Uki of 
both bands and the spin component s|[^ parallel to the 
impurity magnetization (quantization axis) via 




( 20 ) 


As it is common for spin-dependend single particle ener¬ 
gies like the Dresselhaus-®^ or Rashba-terms^®, one could 
expect that the spin-dependend part of the correlation 
energy can be written as an effective magnetic field in 
which the electron spins process. This additional pre¬ 
cession movement could be made responsible for the fre¬ 
quency renormalization discussed above. However, al¬ 
though the corresponding effective field due to the corre¬ 
lation energy has the same form as the frequency renor¬ 
malization, it is larger by a factor of 2. We attribute this 
to the fact that the correlation energy is not an average 
over a Hermitian single particle operator, but comprises 
multiparticle effects, where the naive identification of an 
effective magnetic field can lead to incorrect predictions. 

A particularly interesting and transparent case is that 
where the impurity magentization (S) vanishes. Then, 
the correlation energy takes the form 


{HTdl 


h V 




Am* 

2TTh 


El- 


tdBZ — Wki 
Wki 


nki 


( 21 ) 


Thus, for (S) = 0, we find a logarithmic divergence of 
the correlation energy with respect to ujbz —t oo and 
<^ki 0. In both limits, the correlation energy is neg¬ 
ative and independent of the sign of the coupling con¬ 
stant Jsd- The negative correlation energy suggests a 


formation of correlated carrier-impurity states, i. e., mag¬ 
netic polarons. Such states have also been predicted in 
GaMnAs®®, where the correlation energy stems from an 
attractive Coulomb potential due to the charged state of 
the Mn ions which acts, e. g., as acceptors, if they substi¬ 
tute Ga atoms in the GaAs crystal lattice. An additional 
contribution to the correlation energy originating from 
the spin-dependent s-d interaction could enhance the for¬ 
mation of the magnetic polarons in GaMnAs. Further¬ 
more, our finding is consistent with the fact that also in 
the solution of the Kondo problem, the ground state is 
strongly correlated^. 

Again, the divergence at Wki —>■ 0 in Eq. (21) is in- 
tegrable so that the total correlation energy always as¬ 
sumes finite values. To estimate the magnitude of the 
correlation energy we consider the case where (S) = 0 
and the spectral electron distribution is given by a Gaus¬ 
sian centered at the band edge with standard deviation 
Eg = 1 meV. For the parameters of the Cdo.gsMno.orTe 
quantum well discussed above, we find from Eq. (21) a 
value of w-l.S ^eV per electron. Thus, the cor¬ 

relations can be destroyed by thermal fluctuations when 
the temperature T exceeds 20 mK. 


VI. CONCLUSION 

A microscopic quantum kinetic theory is employed to 
describe the spin dynamics of carriers and magnetic im¬ 
purities in diluted magnetic semiconductors (DMS) ac¬ 
counting also for the dynamics of the carrier-impurity 
correlations. The role of the correlations is examined to 
shed light into the controversy about their importance: 
While some authors assume that the Kondo physics due 
to carrier-impurity correlation is of minor importance®^, 
others®® hnd divergences in a perturbative treatment of 
the spin dynamics in DMS, similarly to the appearance of 
divergences found in the metallic Kondo-effect®. In the 
present study, we hnd that the correlations, besides me¬ 
diating the spin transfer between carriers and impurities, 
are also responsible for a renormalization of the preces¬ 
sion frequency of up to a few percent in two-dimensional 
structures. Since for these values, the relative frequency 
renormalization is negative, the precession frequency of 
the electron spin is reduced. 

In order to experimentally probe the correlation in¬ 
duced frequency renormalization, the spectral features 
of the laser pulse have to be precisely controlled. 
Furthermore, it was reported®® that an antiferromag¬ 
netic impurity-impurity interaction inhuences the ther¬ 
mal equilibrium value of the Mn magnetization, which in 
turn changes the measured electron spin precession fre¬ 
quency. Therefore, it is common to introduce a htting 
parameter Tq and describe the equilibrium Mn magneti¬ 
zation by a Brillouin function with effective temperature 
Teff = To -|- r, where T is the temperature of the sample. 
This complicates the identification of correlation induced 
changes in the precession frequency. To distinguish both 
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effects it is useful that in addition to the dependence 
on the spectral position and shape of the exciting pulse, 
the relative frequency renormalization due to the corre¬ 
lations is independent of the impurity density, while the 
impurity-impurity interaction depends on the mean dis¬ 
tance between the impurity ions and is not influenced 
by the excitation conditions. Because of this and from 
the different parameters entering the prefactor of the fre¬ 
quency renormalization, we find that the most promising 
samples for experimentally accessing the correlation in¬ 
duced frequency renormalization are very narrow quan¬ 
tum wells with large effective masses and a large coupling 
constant Jgd while the impurity concentration should be 
relatively low. Also, we find that the spectral properties 
of ultrashort pulses with durations in the 100 fs range 
suit this purpose. 

Although the k-dependence of the frequency renormal¬ 
ization can in principle lead to a dephasing of carrier 
spins, the spin transfer from the carriers to the impurities 
is usually much faster, so that this dephasing mechanism 
yields only very small corrections to the total decay of 
the carrier spin. 

By comparing the results of the quantum kinetic dy¬ 
namics including the time evolution of correlations with 
the Markov limit, where the correlations can be expressed 
in terms of carrier and impurity variables alone, we find 
that simple Markovian expressions for spin transfer rate 
and frequency renormalization reproduce the results of 
the quantum kinetic treatment very well. Although log¬ 
arithmic divergences are found for some values of carrier 
wave-vectors as well as for bandwidth —>■ oo, similar to 
the metallic Kondo effect, the total frequency renormal¬ 
ization, which is obtained by integrating with a smooth 
spectral carrier-distribution, is always finite. Also, the 
divergences only appear in the Markov limit, while the 
effects induced by the correlations are analytic by con¬ 
struction when a finite memory is accounted for. 

The explicit expressions for the correlations in the 
Markov limit also allow to find an equilibrium value for 


the correlation energy in terms of carrier and impurity 
variables. The form of the correlation energy is simi¬ 
lar to the expression for the frequency renormalization 
and hints towards the appearance of correlated carrier- 
impurity states for low temperatures, independent of the 
sign of the coupling constant. However, the required tem¬ 
peratures are lower than the values that are typically 
considered in experimental setups probing DMS. Thus, 
in these cases, the correlation energy is not expected to 
influence the thermodynamic properties of DMS strongly. 

Although many similarities between the DMS systems 
and the Kondo problem can be found, a one-to-one cor¬ 
respondence can not be established since the density of 
states in semiconductors is much smaller than in met¬ 
als. Hence, the assumption that the number of magnetic 
impurities exceed the number of quasi-free carriers, that 
is valid for DMS and has been used in the derivation 
of the equations of motion on which the present study 
is based^s, is opposite to the situation in usual metallic 
Kondo-systems. Furthermore, the large number of car¬ 
riers in the metal can be expected to lead to a fast car¬ 
rier spin relaxation due to Coulomb scattering. However, 
the approach of finding a quasi-equilibrium value for the 
correlation energy from a quantum kinetic theory could 
be extended and might help to understand the Kondo 
physics from a new point of view. Since we find that the 
correlation effects on the spin dynamics and the correla¬ 
tion energies are linear in the density of states, our the¬ 
ory indeed agrees with the finding of strongly correlated 
states for low temperatures in metallic Kondo systems. 
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